Unifying treatment of nonequilibrium and unstable dynamics of cold bosonic atom 
system with time-dependent order parameter in Thermo Filed Dynamics 
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The coupled equations which describe the temporal evolution of the Bose-Einstein condensed 
system are derived in the framework of nonequilibrium Thermo Field Dynamics. The key element 
is that they are not the naive assemblages of presumable equations, but are the self-consistent ones 
derived by appropriate renormalization conditions. While the order parameter is time-dependent, 
an explicit quasiparticle picture is constructed by a time-dependent expansion. Our formulation is 
valid even for the system with a unstable condensate, and describes the condensate decay caused 
by the Landau instability as well as by the dynamical one. 
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I. INTRODUCTION 

The systems of trapped cold atoms are ideal for study- 
ing the foundations of quantum many-body theories such 
as quantum field theory and thermal field theory. They 
are dilute and weak-interacting, so theoretical calcula- 
tions can be compared with experimental results di- 
rectly. Since the realization of Bose-Einstein conden- 
sates [1-3], many intriguing phenomena have been ob- 
served with good accuracy, and offer opportunities to 
test many aspects of quantum many-body theories in 
both equilibrium and nonequilibrium. Among them, the 
unstable phenomena of the condensate attract our at- 
tention, because firstly to formulate unstable quantum 
many-body systems is still an open problem and and sec- 
ondly nonequilibrium processes accompany the instabil- 
ity in thermal situation. 

Theoretically, the instability of the condensate is char- 
acterized by the eigenvalue of the Bogoliubov-de Gennes 
(BdG) equations [4-6] , which follow from linearization of 
the time-dependent Gross-Pitaevskii (TDGP) equation 
[7]. Since the BdG equations are generally eigenvalue 
ones for non-Hcrmitian matrices , their eigenvalues can 
be complex. The emergence of complex eigenvalues is in- 
terpreted as the sign of the dynamical instability. This 
instability is associated with the decay of the initial con- 
figuration of the condensate and can occur even at zero 
temperature. On the other hand, if the negative eigen- 
values for a positive-norm mode are present, the system 
shows another instability, called the Landau instability. 
This instability, in which the thermal cloud plays an es- 
sential role to drive the condensate toward a lower energy 
state, is suppressed at very low temperature. 

The observations of both the Landau and dynamical 
instabilities are reported in several systems, especially 
in the system where the condensate flows in an optical 
lattice [8, 9], and they are in good agreement with the 
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analyses of the BdG equations [10, 11]. 

Although the TDGP equation outlines the experiments 
corresponding to the dynamical instability at very low 
temperature [13, 14], e.g. the multiply-quantized vor- 
tex splitting [12], an detailed description of the unstable 
dynamics in thermal situations is not trivial. That is be- 
cause there is no more quasi-stable state, and so a fully 
nonequilibrium theory is required. 

There are known two nonequilibrium thermal field the- 
ories, i.e., the closed time path (CTP) formalism [15] 
and the Thermo Field Dynamics (TFD) [16]. The CTP 
formalism is widely used. But we employ the TFD for- 
malism in this paper, because the concept of quasiparti- 
cle picture which is essential for quantum field theory is 
clear even in nonequilibrium situations. In TFD, which 
is a real-time canonical formalism of quantum field the- 
ory, thermal fluctuation is introduced through doubling 
the degrees of freedom, and the mixed state expectation 
in the density matrix formalism is replaced by an average 
of a pure state vacuum, called the thermal vacuum. It 
is crucial in our formulation of TFD to construct the in- 
teraction picture. In quantum field theory, the choice of 
unperturbed Hamiltonian and fields is that of quasipar- 
ticle picture, and concrete calculations are possible only 
when a particular unperturbed representation, or a par- 
ticular particle picture, is specified. One does not know 
an exact unperturbed representation beforehand. 

So far we have investigated the cold atom system 
with a time-independent configuration of the condensate 
in TFD, and derived the quantum transport equation 
which describes the temporal evolution of the quasiparti- 
cle number distribution [19]. It was essential to construct 
an explicit quasiparticle picture there. In contrast to 
the previous investigations [20-24] which are based on a 
phase-space distribution function, our transport equation 
contains an additional collision term which is traced back 
to our choice of an appropriate quasiparticle picture. The 
additional collision term, which we call the triple pro- 
duction term, vanishes in the equilibrium limit if there 
is no Landau instability, but remains non-vanishing to 
prevent the system from equilibrating if there is Landau 
instability. Thus our transport equation with the triple 



production term and the other ones without it predict 
definitely different behaviors of the unstable system. 

In this paper, we derive the coupled equations which 
describe the noncquilibrium dynamics of the cold atom 
system with a time-dependent order parameter. They 
are the TDGP equation, the TDBdG equations, and the 
quantum transport equation. The key points are that 
while the order parameter is time-dependent, we con- 
struct a time-independent quasiparticle picture and so 
that the stable vacuum which are essential for quan- 
tum field theory. These are accomplished by expand- 
ing the field operator with the time-dependent complete 
set evaluating by the TDBdG equations [25] . The quan- 
tum correction to the TDGP equation is determined self- 
consistcntly and simultaneously as the quantum trans- 
port equation by some renormalization conditions [17]. 
Solving the coupled equations numerically, we illustrate 
the dynamics of the condensate decays with either the 
Landau instability or the dynamical one and discrimi- 
nate the two instabilities. 

This paper is organized as follows. We consider the 
cold bosonic atom system with a time-dependent or- 
der parameter at zero temperature in Section II. We 
show that it is crucial to expand the field operator 
by the solutions of TDBdG equations to maintain the 
time-dependent quasiparticle picture. In Section III, 
a nonequilibrium system is considered, and the de- 
grees of freedom are doubled to treat the system with 
the noncquilibrium TFD. We construct a systematical 
method to obtain the coupled equations, and derive those 
explicitly in the leading order. In Section IV, we consider 
a simple system with the Bose-Hubbard model and cal- 
culated the coupled equations numerically. Section V is 
devoted to summary. 



where x — (x,t). Reflecting the existence of the con- 
densate, the field operator ip is divided into a classi- 
cal part C(x) and a quantum one p{x) on the criterion 
(0|(^(a;)|0} = 0. Note that the vacuum is not yet specified 
and that C( x ) is an arbitrary function at this stage, and 
the division must be completed later self-consistcntly. 
The doublet notation is introduced as 
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and the unperturbed Hamiltonian Hq , bilinear and linear 
in ip and ip', is 
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The counter terms 6T(x) and 8C{x) are determined later 
self-consistently. The perturbed Hamiltonian H lnt = H — 
Hq is given as 



II. FORMULATION OF QUANTUM FIELD 
THEORY 

In this section, we consider the cold bosonic atom sys- 
tem with a time-dependent order parameter at zero tem- 
perature. We start with the following Hamiltonian to 
describe the trapped dilute bosonic atoms: 
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where m, V(x), jj,, and g represent the mass of an atom, 
the trap potential, the chemical potential, and the cou- 
pling constant, respectively. Throughout this paper h is 
set to be unity. The bosonic field operator ip(x) obeys 
the canonical commutation relations 
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From the original Heisenberg equation for ip and the 
time-dependent £, the field equation for tp in the interac- 
tion picture should be 

ip = (C + 5T 11 ) p+(M Q + ST 12 ) pt + ho( + 5C-i(. 

(12) 
Due to the last term, this time-evolution is generated not 
by Ho in Eq. (5) but by 



m = Ho-i Id 



& - cv 



as 



vp : 



[<p,K. 



(13) 



(14) 



The condition (0|y(x)|0) = must hold for any t in the 
unpertubed representation, which implies 



i^- t (0\p(x)\0) = (0M0) 



0. 



(15) 



for the time-independent vacuum. According to Eq. (12), 
we have 



5C = {ify - ho)C , 



(16) 



and the time-evolution operator Hq becomes a simple 
quadratic form: 
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A. Field expansion for time-independent order 
parameter 

Before going into further discussions, let us briefly re- 
view the Bogoliubov-de Gennes (BdG) method which 
diagonalizes the unperturbed Hamiltonian Ho (= Hq) 
in case of the time-independent order parameter C,(x). 
The BdG equations are simultaneous eigenvalue equa- 
tions given by [4-6] 



T(x)yt(x) = uJiye(x) . 



(18) 



Since the operator T is non-Hermitian, the eigenvalues 
is not always real but can be complex in general. The 
condition for the emergence of complex eigenvalues in 
the BdG equations has been studied both numerically 
[32-34] and analytically [35-37], and the quantum field 
theoretical formulation has also been discussed [38] . The 
emergence of complex eigenvalues implies the dynamical 
instability of the system, and a drastic temporal change 
of the order parameter occurs then. While our propose 
in this subsection is to find a stable initial condition as 
becomes apparent later, we consider only the case where 
no complex eigenvalue emerges. 

Eigenfunctions belonging to the non-zero real eigenval- 
ues can be orthonormalized under the indefinite metric 
as 



d 3 x y\{x)a 3 ye>{x) = 5 U > , 
d 3 x z\{x)u 3 zt\x) = -bw 
\d 3 xy\(x)o 3 zg,(x) = 0, 



(19) 
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with i-th Pauli matrix <7j. The function ze, defined by 
z t — v\y\-> is an eigenfunction belonging to — ue, when 
yi is an eigenfunction belonging to we- 

Due to the Nanbu-Goldstonc theorem [39], there is a 
zero mode eigenfunction in the BdG equations [30, 31]: 
Tt/o = . The zero mode eigenfunction yo is orthogo- 
nal to all the other eigenfunctions, and what is more, is 
orthogonal to itself. Hence, an additional adjoint mode 
y_i has to be introduced for the completeness as 



where / is determined to satisfy the normalization con- 
dition: 



ld 3 xy' { _ 1 (x)o 3 yo(x) = 1 



(23) 



It is convenient to rewrite the whole orthonormal con- 
ditions with the 2x2 matrix form as 



d 3 x W A (x) W^(x) = 5aa, , 
W t (x) = a 3 f y | (a;) ) a 3 , W~\x) = (y e (x) z e (x)) 
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where 



(25) 
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(26) 

with A = £, 0. 

The completeness condition, 

Yl [yd x )yl( x ') - Zi{x)z\{x') 

+ y (x)y t _ 1 (x r ) + y_ 1 (x)yl(x') = a 3 S(x - x') , (27) 
is simply expressed as 

^W^ix) W A (x') = 5(x - x') , (28) 
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and the field operators in the doublet form are expanded 

as 
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The operators be, p, and q satisfy the canonical commu- 
tation relations [be,b £l ] — Su' and [q,p] = i, respectively, 
and the unperturbed Hamiltonian becomes 



H, 
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Ty_! = Iy , 
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which is diagonalized except for the zero mode sector, 
that is, that of the quantum mechanical operators p and 
q. The choice of the wave function for the sector is not 
settled yet. Although the diagrammatic calculation is 
possible once a particular wave function for the sector 
is chosen, we suppress the zero mode in what follows in 
order to avoid the ambiguity. 



B. Field expansion for time-dependent order 
parameter 

The time-dependent order parameter £(x) implies the 
time-dependent T(x). In order to deal with this situa- 
tion, we suppose a time-dependent orthonormal complete 
set {W^(x)} which is defined by Eq. (25) with some time- 
dependent functions ye{x) and zg(x) and which has the 
properties of Eqs. (24) and (28) at equal time. Then the 
field operators are expanded as 



P a (*) = £w7 1,a/ W(*). 



(35) 



We note that because of the time-dependence of We(x) 
the time-evolution operator for be is not Hq but H$: 



H b (t) = H$(t) 



- jd 3 xd 3 x' (p a {x) 



Y,Wi\x)W i {x') 



a/3 



vV) 



(36) 



Here and hereafter, we take a common time variable for 
x = (x,t) and x' = (x',t). As we are going to treat 
the time-dependent order parameter in the quasi-particle 
picture represented by be, the necessary condition is that 
Hq to be diagonal: 

H b (t) = J2x t (t)bl(t)be(t) (37) 

e 

= 5 £&?(*) [MW^t), (38) 

(. 

where Xe(t) is an arbitrary real function. Therefore, 
We(x) must satisfy 

TWf 1 = iiwf 1 + W/~V 3 . (39) 

We eliminate Xe from the equations by the replacement 

\Yf\x) -> \Vf\x) e^ ds M*)^ and obtain 



T(x)Wf 1 (x) = i^Wf\x), 

or equivalently the TDBdG equations 

d 

T(x)ye(x) = i—y e (x). 



(40) 



(41) 



If the orthonormal complete set {W^(x)} is chosen as the 
initial condition of Eq. (40), it keeps the orthonormal- 
ity and the completeness for all the time because of the 
following equations: 



i— jd 3 x W e (x)W e ' 1 (x) =0, 



'dt 



J2w- 1 (x)W e (x') = Q. 



(42) 
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Thus the procedure of constructing the time- 
dependent complete set is obtained: Solve the eigenvalue 
problem T(x)We(x) — u?eWe(x) at an initial time to, an d 
calculate the time evolution according to i-^We(x) = 
T(x)Wi(x). Next, expand the field operator <p(x) by 
the complete set {W^(a;)} as Eq. (35). This expansion 
is obviously reduced to the ordinary one in the limit of 
time-independent order parameter and has already been 
proposed by Matsumoto and Sakamoto [25]. What we 
have shown in the above paragraph is that the choice of 
We(x) is justified from the viewpoint of quantum field 
theory: the quasi-particle operator be, constructing the 
Fock space on the the Bose-Einstein condensed vacuum, 
diagonalizes the time evolution operator even for the 
time-dependent order parameter. 



III. NONEQUILIBRIUM TFD FORMULATION 

In this section, we double every degree of freedom to 
treat a nonequilibrium system in TFD and introduce the 
thermal Bogoliubov transformation 



with 



,,„, M _ , 1 + ne(t) -ne(t) 



Bf{t) 

i,^/,x _ A ne{t) 
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It is important to take the above particular form of the 
thermal Bogoliubov matrix, as one calls the a = 1 repre- 
sentation [16], which enables us to make use of the Feyn- 
man diagram method in nonequilibrium systems [26] . In 
TFD, the thermal average is represented by the pure 
state expectation of the thermal vacuum, denoted by 
0), and the operators who annihilate |0), are not the 
6-operators but the £-ones: 



&|o> = &|o> = o, <o|eJ = <o|<ej = o 



Ft _ 



(47) 



The number distribution ne(t) is given by 

ne(t) = (0\b\(t)be(t)\0) , (48) 

and its time dependence is determined later. The com- 
bination of the two transformations, £ into b and b into 
ip, involves the 4 x 4-matrix transformations, 



yia _ K -l,/j,ai//3 ^u/3 

Up — D £ Q . 



bf = ^^ , (49) 

^ = j2wr 1 ' l,avP i>? > v* = J2 vt a wF v(3 . ( 5 °) 



with the 4x4 thermal Bogoliubov and BdG matrices 



defined by 



B 



— l,fj,ais{3 



#al£/3l-B 



— 1,(41/ 



+ S a2 Sf32 (o\B t <Ji) , (51) 
BT vfi = 6 al 60iBir + cW,32 {<TiB t ai) liV , (52) 



W, 



— l,(iai//3 



S^W. 



-l,a/3 



(53) 

Wr^ - S*W? , (54) 

where the quartet notations for be are introduced by 



br = 



iy<i 




• ^ 



6J 



6? [ft/ai]" 



- ^t 



(6] -6, -b e b\) 



7t\"P 



(55) 
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and in similar fashions for ^ and ip. The Hamiltonian 
of TFD, which should generate the time translations of 
both the non-tilde and tilde operators, is not the ordinary 
Hamiltonian H but the hat Hamiltonian H = H — H. 
The time-independence of the thermal vacua requires the 
minus sign in front of H. Furthermore, the unperturbed 
Hamiltonian in nonequilibrium TFD is Hq — Hq — Q for 
nonequilibrium system with Hq = Hq — Hq and the ther- 
mal counter term Q(t) , caused by the time-dependence 
oine(t): 



Q(t) 



E^k~™ 



/0 1 



\ 


V 
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if it). 



The unperturbed and full propagators for ip and £ are 



9 il> 1*1) *2 J 



-i(0\T[^ a ( Xl )^(x2)}\0), 
-*<0|T[^(xi)^(x2)]|0), 



v/3. 



-i(0|T[^ Q (il)^7(*: 

^(oir^^K^fe)]^), 



(58) 
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(60) 
(61) 



respectively, where the subscript H denotes a quantity in 
the Heisenberg picture. The self-energies £ and S are 
introduced by the Dyson equations G = A + AUG and 
g = d + dSg, respectively. 

Our critical step is to adopt the following three renor- 
malization conditions simultaneously to determine the 
whole time evolution of the system, explicitly to deter- 
mine the unknown functions ne(t), 5T(x) and £(x) ( or 
6C(x), see Eq. (16)): 



(i) 



(iii) 



3« 21 (M) = 



Re Si 



i— shell 



0. 



\ m (x)\0)=0. 



The condition (i) is what we have proposed for a 
nonequilibrium system with the static condensate [19] as 
a natural extension of the one for non-condensed system 
proposed by Chu and Umezawa [17, 18, 27]. It provides 
the transport equation which determines the temporal 
evolution of the unperturbed number distribution ne(t). 
The possible diagrams in the leading order are indicated 
in Fig. 1. Because the contributions from Fig. 1 (a) and 
(c) vanish, the leading ones come from Fig. 1 (b) and (d), 
and the latter one is proportional to he(t). According to 
the detailed calculations given in Ref. [19], we obtain 



J 



n f 



ht{t) = 4.9 2 Rc / dti^^ {n ei n e2 (l + n e ) - (1 + n £l )(l + ne 2 )nt} ti (ye,Xvv) t (Xvv>Vt) tl 

1 2 + {n ei (l + ni 2 )(l + n e ) - (1 + ne 1 )m 2 ne} ti (ye,Xvz) t {Xvz>ye) tl 

+ {{I + ne 1 )ne 2 {l + ne) - m^l + ne 2 )m} ti (ye,Xzy) t (Xzy,ye) ti 

+ {(l + ne 1 ){l + ne 2 )(l + ne) - ne ± ne 2 ne} ti (ye, Xzz) t (Xzz,ye) ti 

I 



(62) 



where 



x a vy = Cvt y? 2 + Cvt vl + Cvt vt , (63) 

Xy Z = Cvt z? 2 + Cvt z? 2 + Cvt 4 2 > (64) 

x% = Cztvl + Czgyt + Cztvt , (65) 

X a zz = Czt4 2 + Cztzl + Cztzg , (66) 



with C, a (x) 




and a denotes a = 2, 1 for 



a = 1,2, respectively. The subscripts t and t\ are the 
time arguments, and the parenthesis denotes the inner 
product: 



(y,x) 



d 3 xy*' a (x) X a (x). 



(67) 



The term in the fourth line of Eq. (62) is what we call the 
triple production term. As is discussed in Ref. [19], this 
term plays a crucial role if there is the Landau instability, 
but is vanishing due to the energy conservation otherwise. 





(c) 



(d) 



ST 



FIG. 1: One-loop propagators and counter terms. The solid and 
dashed lines denote the unperturbed propagator and the order pa- 
rameter f, respectively. 



While the transport equations with a phase-space distri- 
bution, derived in the other methods previously, lack this 
term, the prediction of our transport equation is different 
from those based on the other transport equations, when 
there is the Landau instability. 

The condition (ii) is the energy renormalization which 
determines ST. Since the leading contribution to the self- 
energy is the tag diagram in Fig. 1(a) and has the form 
of Sw(ti,fa) = S' u ,(ti)5(ti —t 2 ) , the on-shell self-energy 
is naturally defined as S' u {t\) even in the nonequilib- 
rium situation. The concrete forms for the self-energy of 
Fig. 1(a) and (c) are 



S$' Ul \tiM) = S(h -t 2 )(yt,a 3 T^y e ) , (68) 
4 c) ' mi (*i.*2) = S(h-h)(yt,a 3 STye)^ , (69) 



with 



can be rewritten explicitly in terms of ne(t) and yi(x) as 



IvVloHEMlfc 1 ! 



,.2|2\ 



,,2|21 



1 2,* 



W|0) = V(2n£ + l)y^ 



(72) 



(73) 



One can show from Eq. (41) that the time-dependence 
of the total number of non-condensed atoms N cx (t) = 
fd 3 x{0\^(x)tp(x}\0) is 

J t N cx(t) = Y, [( 2n z(t) + l)Im(y e ,T ye) +fn(t)(yt,y t )] 

£ 

(74) 
where 5T(x) drops because of Im{ye,STy^j = 0. Note 
that the second term in Eq. (74) is of two-loop order, 



(b) 



-X ho-idt 



(c) 
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a/3 



T V,afi (x) = a ( 2{0\vHxMx)\0) {0\tp{x)tp(x)\0) , 
1 ' ■" \-(0\<p(xMx)\0)* -2(0\^(x)<p(x)\0)) 

(70) 
Then, we find the condition (ii) is satisfied by 



ST(x) = T^(x) 



(71) 



which is equivalent to the result of the Hartree-Fock- 
Bogoliubov approximation. The matrix elements of T^ 



FIG. 2: Tadpole diagrams. 



although we collect only the one-loop diagrams. This 
comes from the fact that hi is of one-loop order according 
to the quantum transport equation (62), while ng itself 
is of no loop one. 

The last condition (iii) is the self-consistent criterion 
for dividing the original field operator ip into £ and (p. 
The corresponding diagrams at tree and one-loop levels 
are shown in Fig. 2 (a) and (b), respectively. Although 
the energy renormalization is performed at one-loop level, 
the two-loop order correction indicated in Fig. 2 (c) is 
also considered here for the conservations of the total 
number of atoms. The quantity j(x) which is a two- 
loop order part of the counter term SC(x) is determined 
later to cancel the two-loop contribution which appears 
in Eq. (74). The condition with the diagrams in Fig. 2 is 
written as 



d 4 Xl (^ all (x,x 1 ) + ^ a21 (x,x 1 )) 

[h (xi) - id tl + 2.g(0|^ t (x 1 )^( a ; 1 )|0)}c(a:i) + g(0\<p(xi)<p(xi)\0)C(xi) - ij{xi) 
d 4 Xl (A^ Ql2 (a;,xi) + A^ 22 (x,xi)) 

{ho(xi) + id tl + 2g{0\^(x 1 )^{x 1 )\0)}C(x 1 ) + g(0\<p(xi)<p(xi)\0)*C(xi) + 17(31)] = , (75) 



which implies 
i^C = (ho + 2g<0|v?V|0})C + g(0\ w|0)C* - *7 > (76) 



where /iq is defined in Eq. (10). Thus, the modified 



TDGP equation has been derived. To find the function 
j(x), we employ the $ derivative approximation [28, 29], 
which can derive the conserving TDGP equation. In this 
approximation, the self-energy is redefined by the deriva- 
tive of the functional <J> = $[G] as 



E(xi,x 2 ) 



5$ 



SG(xi,x 2 ) 



(77) 



Because the two- loop order modification caused by n^, or 
the thermal counter term Q in other words, is required 
here, we only calculate the part of $ which is related to 
Q in the leading order. The contribution of the thermal 
counter term to the self-energy is 



S^^,^)^^^ 1 ^!)^ 1 ^) 



x S QA e 2 {t 1 ,t2)Be 2 {t2)yVe 2 (x 2 ) 



fiavft 



with 



iau/3 



£q^! 2 (*i>*2) = -int 1 (t 1 )S(t 1 -t 2 )5i 1 i 2 
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and that to the functional $ is 



$ Q = d 4 Xl d 4 x 2 ^ f3 (x ll x 2 )A-^ a (x 2 ,x 1 )e'"e a 



^2 dtnt(t)(ye,ye) t , 



(80) 

(81) 



with the sign factor, e 1 — 1 and e 2 = — 1. According to 
the $ derivative approximation, j(x) which is the two- 
loop correction involved in Sq is found to be 



l(x) = 



5<f>Q 
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(79) By substituting Eq. (62) into this, we obtain 



(82) 



7 (x) = ,g 2 Re 



dt i 



E 



(yt 3 ,ye 3 



{n in2 (l + n 3 ) - (1 + ni )(l + n 2 )n 3 } ti Xl yz (x)(xyy,y 3 ) ti 
+ {m(l + n 2 )(l + n 3 ) - (1 + ni)n 2 n 3 } ti X 4 zz (x)(x V z, y 3 ) ti 
+ K 1 + ni)n 2 (l + n 3 ) - m(l + n 2 )n 3 } ti Xl yz (x)(xz V , y 3 ) ti 

+ K 1 + ni)(l + n 2 )(l + n 3 ) - niften 3 } tl X^.(x)(x**, ys) ti 



r 



(83) 



where 



X" 



ft ft ft i ft ft ft i ft ft ft 

2/f J/2 4 + ViV^s + ViV 2 z 3 , 



(84) 



and in similar fashions for X" z , X° '■ and X^ zz . Then, 
the time derivative of the total number of condensed 
atoms No(t) — jd 3 x \((x)\ 2 becomes 

— N a = -^ [{2ni + l)Im(y£,T y £ ) + ht{y i ,y e )] , 

which cancels Eq. (74) and implies the conservation of 
the total atom number. 

Thus, we obtain the coupled equations which describe 
the temporal evolution of the condensed system, those 



are the TDGP equation (76), the TDBdG equations (41), 
and the quantum transport equation (62). 



IV. NUMERICAL RESULT 

In this section, we calculate the coupled equations nu- 
merically to illustrate the nonequilibrium dynamics, es- 
pecially the condensate decays with either the Landau 
instability or the dynamical one, and confirm the quali- 
tative difference between both the instabilities. For this 
propose, we consider a very simple system with the one- 



dimensional Bose-Hubbard model with the Hamiltonian 



H = 



/^ 



U 



-Jlpj {lp l+1 + Ipi-i} - fllpjlpi + -7rlpjlpj4>ilp 



( 86 ) 
Here J, U, and i represent the inter-site hopping, the on- 
site couping, and the site index, and we put the number 
of sites I s — 21, the total number of atoms N = 210, and 
U/J = 0.05. The condensate is introduced as ipi = Q+fi 
with the criterion (0|</3i|0) = 0. It is straightforward to 
apply the method developed in the previous section to 
this model and to derive the coupled equations for this 
system in the leading order. 

To illustrate the condensate decay with the Landau 
instability or the dynamical one, we consider the follow- 
ing situation. First, the equilibrium state with no con- 
densate flow at a temperature To is prepared. Then, 
the condensate is forced to flow instantaneously with the 
quasimomentum k : Q — > Qe lkXi , and the system turns 
into nonequilibrium. This nonequilibrium state is cho- 
sen as the initial state of the calculation, and then the 
coupled equations are calculated numerically. 

Solving the BdG eigenvalue equations at zero tempera- 
ture analytically, we obtain the stability diagram as Fig. 3 
and find that the system is stable for kL/2ir = and 1, 
Landau unstable for kL/2ir = 2 to 5, and dynamically 
unstable for kL/2ir = 6 to 9. Although the result is for 
the zero temperature, it is also expected to be valid for 
the nonequilibrium case with a sufficiently small initial 
depletion. 




0123456789 



kL/2-K 

FIG. 3: Stability diagram for U/J = 0.05 obtained analytically 
from BdG eigenvalue equation at tree-level. The symbols k and q 
stand for the quasi-momenta of the condensate and the excitation, 
respectively, and L is the length of the system. The white, red, and 
black colored cells denote stable, Landau unstable, and dynamically 
unstable regions, respectively. 

The coupled equations are calculated numerically with 
several initial conditions. The depletion of the conden- 
sate is indicated in Fig. 4, and the stable and unstable 
behaviors are clearly discriminated. For the stable condi- 
tion kL/2ir — 1, the nonequilibrium depletion oscillates 
around the equilibrium one and never grows. The deple- 
tion grows initially with oscillation for the Landau unsta- 
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FIG. 4: 



Temporal evolution of the depiction -^ 2Zi (0|vl Vi |0). 
Solid (dashed) line indicates the result with smaller (larger) ini- 
tial depletion, corresponding to a lower (higher) To. Blue, red, and 
black lines correspond to the stable, the Landau unstable, and the 
dynamically unstable conditions, respectively. 
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FIG. 5: Lifetime of the condensate vs initial depletion. Red and 
black lines correspond to Landau and dynamically unstable condi- 
tions, respectively. The lifetime is simply denned by the time when 
the depletion reaches to 0.5. The nonmonotonicity of the lifetime 
in case of the Landau instability is due to the oscillating growth of 
the depiction. 



ble condition kL/2ir = 5, but grows much more rapidly 
and exponentially from the beginning for the dynamically 
unstable one kL/2ir — 6. The decay speed with the Lan- 
dau instability tends to increase if the initial depletion 
becomes large as is shown in Fig. 5, reflecting the fact 
that the Landau instability is caused by the collision be- 
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FIG. 6: The effect of the triple production term to the depletion 
growth. Solid line is the same as in Fig. 4, while doted line indicates 
the result, obtained when the triple production term is omitted. 



tween condensate and non-condensate particles and that 
the collision becomes more frequent for larger depletion. 
On the other hand, the value of the initial depletion is 
not relevant for the dynamical instability, since the non- 
condcnsatc particle plays no essential role then. 

The effect of the triple production term which is a dis- 
tinguishing feature of our quantum transport equation 
is shown in Fig. 6. No qualitative difference between 
the cases with the triple production term and without it 
is notable in both the stable and dynamically unstable 
conditions. That is because that the triple production is 
suppressed due to the energy conservation in the former, 
and that the rapidly growth is governed by the TDGP 
and TDBdG equations but not by the transport equation 
in the latter. In the Landau instability condition, on the 
other hand, the difference is remarkable. Basically, the 
growth behavior disappears if the triple production term 
is omitted. Thus we conclude the essential contribution 
of the triple production term in describing the Landau 
instability. 

Although we present only the results of the three typi- 
cal values of k, representing the stable, Landau unsta- 
ble, and dynamically unstable situations, respectively, 
the qualitatively similar results are obtained for any al- 
lowed value of k. 



V. SUMMARY 

In this paper, the self-consistent equations which de- 
scribe the nonequilibrium dynamics have been derived 
by applying the nonequilibrium TFD to the condensed 
atom system. The system with a time-dependent order 
parameter is considered as the further extension of our 
previous study where the stationary order parameter has 
been assumed. We only have been able to describe the 
first stage of the condensate decay with the Landau in- 
stability heretofore. Now we can predict the next stage 
of the decay dynamics, and can describe not only the 
Landau instability but also the dynamical one. 

To treat the time-dependence of the order parameter 
within the time-independent quasiparticle picture, the 
field operator is expanded with the complete set evalu- 
ating by the TDBdG equations. This method, proposed 
first by Matsumoto and Sakamoto without a detailed vin- 
dication, is obtained here to construct the particle oper- 
ator time- independent. The renormalization conditions, 
the the Chu-Umezawa's diagonalization condition (i), the 



energy renormalization (ii), and the criterion for diving 
the non-condensate and the condensate (iii), are applied 
to determine the coupled equations which are the quan- 
tum transport equation, the TDBdG equations, and the 
TDGP equation. 

The point is that the coupled equations we obtained 
are not the naive assemblages of presumable equations, 
but are the self-consistent ones derived by the appro- 
priate renormalization conditions. To confirm that our 
coupled equations can describe both Landau and dynam- 
ical instability, we consider an one-dimensional Bosc 
Hubbard model, and calculate numerically the depletion 
whose growth implies the decay of the condensate. Since 
the initial depletion is sufficiently small, the unstable 
condition is well characterized by the eigenvalues of the 
BdG equations at tree-level. We found that the depletion 
growths with both the Landau and dynamically unstable 
condition, while dose not with stable condition. Pre- 
dictably in the Landau unstable condition, the growth 
is slower and has stronger initial depletion dependence 
than that in the dynamically unstable condition. 

As we reported previously, our transport equation con- 
tains the triple production term which is absent in the 
one of the other methods. The difference originates in 
the choice of the quasiparticle picture which is essential 
for the quantum field theory. We construct the quasi- 
particle faithfully to the quantum field theory, while the 
others are based on a phase-space distribution function 
and no explicit particle representation is given. This dif- 
ference is inconspicuous if there is no Landau instability, 
but causes a great qualitative change if there is Landau 
instability. Basically, we find that the condensate dose 
not decay without the triple production term in the Lan- 
dau unstable condition, although the decay has been ob- 
served experimentally. We emphasize that the choice of 
the quasiparticle picture is even more essential for a un- 
stable case. A quantitative comparison to a experiment 
of the condensate decay will be the future task. 
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